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This note focuses on the development of a 2D model of a thin liquid film flowing down inside a 
vertical pipe. This model is based on the large wavelength assumption and valid for high Reynolds 
and Weber numbers. 

Un modele 2D decrivant le ruissellement d'un film liquide mince a I'interieur d'une conduite 
cylindrique maintenue verticale est developpe. II est base sur I'hypothese de grande longueur d'onde 
et valable pour des grands nombres de Reynolds et de Weber. 

I. INTRODUCTION 

Various models relate to falling films on vertical or inclined plate and small Reynolds number whereas most exper- 
imental observations at high Reynolds number are made inside cylindrical tubes. The aim of this note is to present 
a 2D model which accounts for small and medium cylinder curvature at high Reynolds number, in order to make a 
more realistic comparison with the experimental data 

II. MAIN EQUATIONS 

We consider a viscous fluid film flowing down a vertical infinite pipe of radius Rc having a common interface 
with a quiescent gas phase. System 111 1411 defines the full equations of motion for the annular falling film. Two small 
parameters are defined : e which is the ratio between the film thickness /iq and the characteristic length and for 
the reduced curvature of the cylinder. In order to reduce these model equations, we use a boundary-layer theory. To 
satisfy the assumptions of large Reynolds numbers (ei?e = 0(1) ) and e^W = 0(1), equations Hlfllll|l are simplified 
by dropping the terms of order e/ Re or higher. We assume that the instantaneous streamwise velocity u is described 
by a self-similar velocity profile Ijl4|l . These reduced equations are integrated through the local and instantaneous 
thickness h of the Hquid film. The final model is a system of two coupled equations iflfih . ITfih l involving h and the 
local fiow rate g, as the only dependent variables. We give an approximated model at order e^. 

III. DISCUSSION 

If one follows the Shkadov procedure [i| the smallness parameter e may be effectively included in the dimensionless 
numbers of the fiow, avoiding to specify a characteristic length Ic- The fact remains that, in the cylindrical geometry 
the parameter associated with the curvature cannot be eliminated in the same way and the specification of Ic is 
mandatory. In this case, the characteristic length may be defined as the length of the neutral mode or the length 
of the most energetic wave P|. The cylindrical model developed in this study could be improved by writing the 
streamwise velocity field as polynomials functions Q. The dimensionless equation are found by using time and space 
scales based only on the physical properties of Hquid which are the kinematic viscosity and the gravity acceleration. 
This may contribute to give more information of the nonlinear terms included in the functions $1 and $2- 
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I. INTRODUCTION 



Les divers developpements theoriques realises sur le film tombant se divisent en deux groupes. Le premier prend en 
compte les deux phases de I'ecoulement a travers des equations d'evolution decrivant les variables de chaque phase. Le 
second groupe, dans lequel nous nous plagons, considere que la phase gazeuse est inactive. Cependant, les difi^erents 
modeles Q existant dans la litterature ne concernent que le film ruisselant sur une plaque plane verticale ou inclinee 
alors que la majeure partie des observations experimentales 0, compatibles avec des situations industrielles, sont 
faites a I'interieur des tubes cylindriques de diametre variant de 0.5 a 10 cm. Pourtant, on salt par experience que 
le cylindre est plus stable a haut debit que le cas plan ce qui n'est pas explicable par la seule periodicite azimutale. 
L'objet de cette note est de presenter un modele 2D qui tienne compte de la courbure du cylindre et les grands 
nombres de Reynolds afin de permettre une future comparaison plus reahste avec les donnees experimentales Q. 



II. PRESENTATION DU PROBLEME 



Nous considerons le ruissellement d'un film liquide mince le long de la parol interne d'une conduite cylindrique 
verticale (figurenj de section circulaire et de rayon Rc sous I'effet de la gravite. Le Hquide est newtonien incompressible 
et isotherme. Les points de depart de la recherche des equations d'evolution du film sont les equations de Navier-Stokes : 
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Fig. 1 - Representation de la geometrie de I'ecoulement. 



V ■ u= 0, (la) 



d u 



p y^Qf+ ■ u j = p g - Vp + A^Aw, (lb) 

ou p est la densite, ^ la viscosite dynamique, p la pression et u= (w, u) la vitesse. A ces equations sont jointes les 
conditions aux limites : 

- de non glissement a la parol : 

u = w = 0, (2) 

- cinematique a I'interface : 



(3) 



3 



- de saut dynamique a la surface libre 

P-Pg = 



-1 3/2 



Rc 
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ou a est le coefficient de tension superficielle et Pg la pression du gaz supposee constante. 

La solution de base correspond a un film tombant avec une epaisseur constante ho le long de toute la conduite. 

d 

Cette solution est obtenue a partir des equations precedentes en supposant que I'ecoulement est permanent ('q^ ~ ^) 

d 

et etabli suivant I'axe (v — 0, — — 0) de la conduite, dans ce cas nous obtenons des equations reduites dont la 

dx 

solution est la composante axiale de la vitesse 
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ou v est la viscosite cinematique. Le debit est 



correspondant a une vitesse debitante 
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L'epaisseur ho est solution de I'equation 



8 Re {rI - [Rc - hof^ = g ho R^x 
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1 /3 

ou Re — uo ho/v est le nombre de Reynolds. L'ecart entre ho et l'epaisseur de Niisselt (3 Re/g) est negligeable 
pour Rc > 1.5 cm. 

Dans la suite, nous introduisons une nouvelle variable y = Rc — r, I'expression © de la vitesse devient 
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correspond au profil de base d'un film liquide ruisselant sur un plan vertical et 
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1 — — - I log 1 — — - I est le terme de correction du profil plan lie a la courbure du cylindre. La parol du 
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2v 

cylindre est identifiee par I'equation y = et I'interface par la relation y — h{x,t) 

Pour tenir compte de I'hypothese de grande longueur d'onde ^.4], nous introduisons le rapport d'aspect e = ho/ Ic ou 
Ic est la longueur caracteristique suivant I'axe des x. Nous negligeons les termes d'ordre superieur a pour aboutir 
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aux equations reduites de Navier-Stokes sous la forme adimensionnelle : 
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La condition d'adherence et les conditions aux limites restent inchangees alors que les conditions dynamiques donnent 
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ou W — a/ (p/iqUq) est le nombre de Weber, R = Rc/h^ le rayon adimensionnel et Fr — uo/y/gho le nombre de 
Froude qui prend la valeur y/Re/3 en geometrie plane verticale. Dans les equations ltTol- [TT|l . les differentes variables 
ont ete reduites de la maniere suivante : 



X -> x/lc, r ^ R- y/hf), t t/{lc/u„), u u/uq, v ev/u^, p p/{pul). 



(12) 
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II est a noter que -77 ^ ^ e. L'hypothese de grande longueur d'onde impose e petit devant I'unite. En faisant 

dt dx 

l'hypothese que eRe — 0{1) et e^W = 0(1) , nous pouvons negliger dans les equations IjlOb . llOb. Illh. ITTbl le terme 
en e/Re, / Re et e^, pour obtenir : 



du dv 



dx dy R~ y 







1 du d^u 
Re \ r dy dy"^ 



1 



du du du dp^. p 
dt dx dy dx ^ 



dp 
dy 



= 



P\y=h =Pg 



W 

R-h 



d^h 



dx^ 2{R-h) \dx J 



dhV 



du 
dy 



\y=h = 0. 



(13a) 
(13b) 
(13c) 

(13d) 

(13e) 



III. EQUATIONS D'EVOLUTION DU FILM TOMBANT 

II apparait done naturel de ramener toutes les equations reduites ifTTH -ITTll a I'interface par integration sur I'epaisseur. 
Cette technique est d'usage courant dans les problemes ou I'inconnue principale est localisee sur une partie de la 
frontiere du domaine. Le relevement des equations IjlOb . llOb l sur I'epaisseur se fait avec l'hypothese que le profil de 
Vitesse est auto similaire 

u{x,y,t) = U{x,t) ^[^^ (14) 
ip{h{x,t}} 

ou tp{y) = y/R— 1/2 {y/R)^ + (1 ^ h/R)^ In (1 — y/R) et U la vitesse locale a I'interface du film. On pose i^iy) = 
(p{y) /(p {h {x, t)). A partir de la definition de q — J ^1 — u dy, le debit instantane du film, le relevement integral 
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consiste a ecrire ^ll^h . \ll^b ) sous la forme 
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En tenant compte des conditions aux limites H13I11II . on aboutit done a un systeme de deux equations dont les 
inconnues sont h I'epaisseur instantanee et q : 
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f'^ du 
ou / (R — y) u^dy = Rq^^Ah), — — L=o = Q^2{h)- Les expressions des fonctions $1 et ^2 sont facilement 
Jo dy 



determinees par les relations 



$1 (h) = R' 



(C) CdC]/[ ^ (C) CdC et $2 (h) = R^'iO)/ I ^ (C) CdC 



En definissant un nouveau petit parametre e,. = qui mesure les effets de courbure, nous faisons un developpement 
en serie des equations lfTfit i-[T6h') et obtenons : 
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--Merh) = i + V0^r+'^^h4, 5M) = -^-|^^-|^^, Cie..h)^2h-hh^. 

On retrouve le cas plan en faisant tendre vers I'infini le rayon de la conduite, ce qui redonne le modele de Shkadov 
correspondant aux equations ifTTh lTTh) avec — 0. 



IV. DISCUSSION ET CONCLUSION 



Ce modele decrivant le ruissellement d'un film liquide dans une conduite cylindrique verticale permet d'explorer 
la dynamique lineaire du film dans des gammes de nombres de Reynolds correspondant a un regime d'ecoulement 
interne laminaire (^ 700). 

Si on suit la demarche de Skhadov P| on pent certes inclure le parametre de petitesse dans les nombres adimen- 
sionnels de I'ecoulement et ainsi eviter de specifier une echelle de longueur axiale. Ce n'est plus le cas en geometric 
cylindrique ou les termes de courbure imposent d'expliciter Ic- Dans ce cas, on peut definir la longueur caracteris- 
tique comme etant la longueur du mode neutre Q ou la longueur d'onde de I'onde la plus energetique Une autre 
demarche consiste d'une part a prendre comme echelle de longueur caracteristique Ic construite a partir de deux 
grandeurs physiques telles que et 5 Q et d'autre part a introduire un nouveau petit parametre lie a la courbure du 
cylindre dont I'expression sera Ic/Rc- Cette demarche permettrait d'obtenir un modele dont la construction se ferait a 
partir du profil de vitesse semi-parabolique dont la forme polynomiale s'adapte bien aux developpements deja inities 
par Q et d'avoir acces aux termes en et ceux propres a la configuration cylindrique et qui sont englobes dans $1 et 
4>2. En pratique le film plan est repute moins resistant aux perturbations exterieures (vibrations, action de I'ecoule- 
ment gazeux,...) que le film annulaire. Pourtant, du point de vue theorique, en considerant la tres faible infiuence du 
rayon de courbure sur I'epaisseur moyenne du film, il est generalement considere que son role est negligeable. II serait 
done interessant, en partant d'une analyse lineaire des instabilites axiales, de quantifier I'infiuence du rayon sur les 
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taux de croissance. 
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